JOURNAL OF GEOPHYSICAL RESEARCH: ATMOSPHERES, VOL. 118, 6505-6520, doi:10.1002/jgrd.50313, 2013 


Impact of surface roughness and soil texture on mineral dust emission 
fluxes modeling 

Laurent Menut, 1 2 3 4 5 Carlos Perez, 2,3 Karsten Haustein, 4,5 Bertrand Bessagnet, 6 
Catherine Prigent, 7 and Stephane Alfaro 8 

Received 25 October 2012; revised 27 February 2013; accepted 2 March 2013; published 22 June 2013. 

[i] Dust production models (DPM) used to estimate vertical fluxes of mineral dust 
aerosols over arid regions need accurate data on soil and surface properties. The 
Laboratoire Inter-Universitaire des Systemes Atmospheriques (LISA) data set was 
developed for Northern Africa, the Middle East, and East Asia. This regional data set was 
built through dedicated field campaigns and include, among others, the aerodynamic 
roughness length, the smooth roughness length of the erodible fraction of the surface, and 
the dry (undisturbed) soil size distribution. Recently, satellite -derived roughness length 
and high-resolution soil texture data sets at the global scale have emerged and provide the 
opportunity for the use of advanced schemes in global models. This paper analyzes the 
behavior of the ERS satellite-derived global roughness length and the State Soil 
Geographic data base-Food and Agriculture Organization of the United Nations 
(STATSGO-FAO) soil texture data set (based on wet techniques) using an advanced DPM 
in comparison to the LISA data set over Northern Africa and the Middle East. We explore 
the sensitivity of the drag partition scheme (a critical component of the DPM) and of the 
dust vertical fluxes (intensity and spatial patterns) to the roughness length and soil texture 
data sets. We also compare the use of the drag partition scheme to a widely used 
preferential source approach in global models. Idealized experiments with prescribed 
wind speeds show that the ERS and STATSGO-FAO data sets provide realistic spatial 
patterns of dust emission and friction velocity thresholds in the region. Finally, we 
evaluate a dust transport model for the period of March to July 2011 with observed 
aerosol optical depths from Aerosol Robotic Network sites. Results show that ERS and 
STATSGO-FAO provide realistic simulations in the region. 
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1. Introduction 

[ 2 ] Dust aerosol particles are produced by wind erosion in 
arid and semiarid surfaces. From the local to the global scale, 
modeling of mineral dust remains a challenge. While dust 
distribution and dust effects are important at global scales, 
they strongly depend on dust emission, which is a threshold, 
sporadic and spatially heterogeneous phenomenon, locally 
controlled on small spatial and temporal scales. Since dust 
plays a key role in regional to global climate and air quality 
because of its effects upon radiation, ocean biogeochemistry, 
and human health, it is necessary to accurately represent its 
emission, transport, and deposition at different spatial and 
temporal scales. 

[ 3 ] During the last years, many efforts have been devoted 
to the reduction of emission uncertainties. Current theo- 
retical knowledge would allow a satisfactory calculation 
of the vertical dust flux in models if the required input 
parameters — surface, soil, and meteorological features — 
were accurately determined. However, the application of 
complex emission schemes in global and to a lesser extent 
regional models is mainly hampered by the lack or the strong 
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uncertainties of the required input data at the pertinent scales 
and the inaccuracies of the driving meteorological model 
mainly with respect to surface wind velocities, turbulence, 
and stability. 

[ 4 ] In order to partly overcome these limitations, current 
global models [e.g., Tegen et al., 2002; Mahowald et al., 
2003; Miller et al, 2006] have assumed varying degrees 
of simplification in the dust emission scheme as a function 
of the availability and accuracy of the input data, and in 
most cases several parameters are tuned to match quanti- 
tative dust observations that are mainly available far away 
from sources. 

[ 5 ] Regional mineral dust models have been also 
developed and used for particular regions such as, for exam- 
ple, DREAM [Nickovic et al., 2001; Perez et al, 2006] 
and CHIMERE [Menut et al, 2007] in Northern Africa. 
Although simplifications also apply to regional models, 
some of them include specific surface and soil data sets 
developed for the most well-known source regions [e.g., 
Chatenet et al, 1996; Callot et al, 2000]. Numerous sensi- 
tivity experiments have been performed to better understand 
and quantify strengths and weaknesses of different formula- 
tions and data sets. The sensitivity to various dust production 
models and meteorology was quantified in Todd et al. [2008] 
over the Bodele region and Menut [2008], Cheng et al. 
[2008], and Zhao et al. [2010], among others. 

[6] In dust model intercomparison exercices [Textor et al , 
2006, 2007; Todd et al., 2008; Huneeus et al, 2011] surface, 
soil, meteorology, and dust production models (DPMs) vary 
from one model configuration to the other. It is therefore dif- 
ficult to attribute model differences to a specific parameter 
or forcing, and error compensations are always possible. In 
an interesting effort, Darmenova et al. [2009] made an in- 
depth comparative analysis of two different physically based 
schemes which were originally developed by Marticorena 
and Bergametti [1995] and Shao et al. [1996]. The relative 
importance of the input parameters was assessed, and some 
recommendations were provided on the selection of input 
parameters, including land and meteorological variables, to 
achieve an improved modeling of dust emission in Cen- 
tral and East Asia. A more recent study implemented three 
dust schemes in WRF-Chem and compared the vertical dust 
fluxes, highlighting large differences for the same meteo- 
rological conditions and surface and soil data [ Kang et al., 
2011]. In parallel, Wang et al. [2012] implemented and com- 
pared two emissions schemes in Community Multiscale Air 
quality Modeling System (CMAQ). 

[ 7 ] In this paper, the sensitivity of mineral dust fluxes to 
surface and soil features is studied over Northern Africa. 

[8] Concerning the soil, the fractions of clay, silt, and sand 
from global soil data sets [e.g., Zobler, 1986] are generally 
mapped according to the well-known textural triangle which 
is based on measurements using wet sedimentation tech- 
niques that may strongly overestimate the amount of loose 
clay particles [. Laurent et al., 2008]. In this paper, we will 
compare this approach against an alternative approach using 
dry techniques [ Chatenet et al, 1996]. 

[ 9 ] Another very important parameter is the aerodynamic 
roughness length, z 0 , that controls the threshold wind fric- 
tion velocity above which saltation starts. The estimates 
of roughness length used in meteorological models do not 
apply to the scales needed for describing the dust emission 


process. Only recently, new encouraging approaches based 
on satellites [ Prigent et al, 2005; Laurent et al, 2008] were 
proposed for providing the information on aerodynamic 
roughness lengths appropriate for global and regional dust 
models. We test and compare the spatial distribution of emis- 
sions using a drag partition scheme with different roughness 
length estimates [ Prigent et al, 2005; Marticorena and 
Bergametti, 1995]. 

[ 10 ] Partly due to the lack of aerodynamic roughness 
length data sets at a global scale, a very common approach 
has been the prescription of an erodibility factor that 
accounts for the spatial distribution of dust source inten- 
sities. Prospero et al. [2002] showed that enclosed basins 
containing former lake beds or riverine sediment deposits 
are preferential sources dominating global dust emission. 
Several model representations of preferential sources have 
been used based on topographic [Ginoux et al, 2001], 
hydrological [Tegen et al, 2002], or geomorphological 
[Zender et al, 2003] approximations. Also, surface 
reflectance retrieved from Moderate Resolution Imaging 
Spectroradiometer (MODIS) [Grini et al, 2005], the fre- 
quency of high Aerosol Index values [ Westphal et al, 2009], 
and the UV- visible surface albedo [Morcrette et al, 2009] 
have been used to identify preferential sources. The results 
of sensitivity experiments with different preferential source 
fonnulations at a global scale [Zender et al, 2003; Cakmur 
et al, 2006] are rather inconclusive since they are depen- 
dent on the model and evaluation data as well as on the 
optimization and evaluation methods. The most popular rep- 
resentation of sources in dust modeling which has shown 
significant improvements is the topographic approach of 
Ginoux et al. [2001]. In this paper, we also test and ana- 
lyze this approach in comparison to the use of aerodynamic 
roughness lengths in a drag partition scheme. 

[ 11 ] The paper is structured as follows. Section 2 describes 
the surface and soil data sets compared in this study. 
Section 3 describes the Dust Production Model (DPM) 
used. The sensitivity studies are performed using a combi- 
nation of idealized and real-world model experiments: In 
section 4, we quantify and discuss, under idealized mete- 
orology, the spatial distribution of emissions in the region 
when using different formulations and data sets describing 
the soil and surface conditions in the region, and in section 5, 
we use a dust transport model to compare real-world sim- 
ulations with different formulations against aerosol optical 
depth measurements in the region. A summary is given 
in section 6. 

2. Surface and Soil Data 

2.1. Soil Type and Soil Texture Data 

[ 12 ] In this study, we use two different soil data sets, 
hereafter named LISA and STATSGO-FAO. The regional 
LISA-soil data set was developed by [Marticorena and 
Beigametti, 1995] for Northern Africa. The data are 
available at http://www.lisa.univ-parisl2.fr/mod/data/. The 
STATSGO-FAO is a global data set used in meteorological 
modeling and which may be used for global and regional 
dust transport modeling, and the data are available at http:// 
soils.usda.gov/survey/geography/. 

[ 13 ] Both provide soil texture characteristics over 
longitude/latitude regular grids. The native resolution of 
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Figure 1. The LISA and STATSGO-FAO soil types data 
reggridded to 1° x 1° resolution. The values correspond to 
the surface types explained in Table 1. For each cell, only 
the most important relative surface is displayed (on a total of 
five possible surface types in a cell). 


the data sets are different: 1° x 1° for LISA soil and 
0.0083° x 0.0083° for STATSGO-FAO. In order to have 
homogeneous data sets, the STATSGO-FAO data are regrid- 
ded into a regular 1° x 1° grid, retaining the dominant 
feature as presented in Figure 1. 

[ 14 ] The soil types in the two data sets are listed in Table 1 
together with their relative percentage of clay, silt, salt, and 
sand (coarse and fine-medium). Note that STATSGO-FAO 
does not contain salts and that fine-medium sand and coarse 
sand are distinguished according to criteria detailed in Tegen 
et al. [2002]. With this approach, clay loams are considered 
as highly unlikely to contain coarse sand, and sandy clay 
loam can contain both coarse and medium fine sand. The 
textural triangle is based on measurements performed by wet 
sedimentation techniques which break the soil aggregates, 
leading to high amounts of loose clay particles that generally 
form aggregates of larger size and that may not be encoun- 
tered in natural soils. Following Tegen et al. [2002], clay is 
in the form of aggregate in loamy sands, and it is reassigned 
to the silt fraction. In the case of LISA, only the clayey soil 
type contains clay particles. 

[ 15 ] Each soil population can be described in tenns of 
a mean mass median diameter, D p (p m), and an associated 
standard deviation a (Table 2). An issue is that the two 
data sets do not use the same values. For the LISA data 
set, the characteristics of the soil types are from Chatenet 
et al. [1996], and for the STATSGO-FAO from Tegen et al. 
[2002]. In this study, for the sake of comparison, we use the 
same size -distribution characteristics for each population in 
the soil as displayed in the third column of Table 2. 

[ie] Each LISA-soil grid cell contains the proportion of 
the five most dominant soil types. For the highly resolved 
STATSGO-FAO data set, we follow the same approach by 
counting the relative proportion of the five most dominant 
soil types in 1° x 1° cells. 


Table 1. Characteristics of Each Model Soil Type 3 


Number CS 

FMS 

% 

Silt 

Clay 

Salt 

Soil Type 



LISA data soil types 


i 

0 

62.5 

37.5 

0 

0 

silty fine sand 

2 

80 

20 

0 

0 

0 

medium sand 

3 

100 

0 

0 

0 

0 

coarse sand 

4 

90 

10 

0 

0 

0 

coarse medium sand 

5 

0 

100 

0 

0 

0 

fine sand 

6 

31.25 

31.25 

37.5 

0 

0 

silty medium sand 

7 

0 

0 

20 

0 

80 

salty medium sand 

8 

0 

0 

8 

0 

92 

salty silt fort 

9 

0 

0 

0 

100 

0 

clayey soil 

10 

0 

0 

50 

0 

50 

salt waste 

11 

0 

0 

100 

0 

0 

agricultural soil 

12 

0 

50 

10 

0 

40 

salty fine sand 

13 

60 

0 

40 

0 

0 

silty coarse sand 



STATSGO-FAO data soil types 

1 

46 

46 

5 

3 

0 

sand 

2 

41 

41 

18 

0 

0 

loamy sand 

3 

29 

29 

32 

10 

0 

sandy loam 

4 

0 

17 

70 

13 

0 

silt loam 

5 

0 

10 

85 

5 

0 

silt 

6 

0 

43 

39 

18 

0 

loam 

7 

29 

29 

15 

27 

0 

sandy clay loam 

8 

0 

10 

56 

34 

0 

silty clay loam 

9 

0 

32 

34 

34 

0 

clay loam 

10 

0 

52 

6 

42 

0 

sandy clay 

11 

0 

6 

47 

47 

0 

silty clay 

12 

0 

22 

20 

58 

0 

clay 


a CS = coarse sand; FMS = fine-medium sand. 


2.2. Land Use Data 

[ 1 7 ] Land use data describe the surface occupation (wood, 
mountains, lakes, etc.) and is typically used to prescribe 
some properties in meteorological models, such as the 
roughness length, z 0 . For each land use type, a z 0 is typi- 
cally prescribed. As described in the next section, roughness 
lengths based on land use types are not appropriate for use 
in dust emission parameterizations. 

[is] For mineral dust emission calculations, the land use 
is typically used to provide a desert mask specifying what 
surface is potentially erodible. The LISA data set already 
contains this information as the percentage of erodibility 
in each grid cell. Global or regional models typically con- 
sider arid and semiarid categories from a land use database 
to identify potential dust source areas. Here, to comple- 
ment the STATSGO-FAO soil data set, we use the U.S. 
Geological Survey (USGS) land use data set with a high hor- 
izontal resolution (0.0083° x 0.0083°). As in the case of the 
STATSGO-FAO data set, we calculated the relative surface 

Table 2. Mean Characteristics of the Soil Textures 3 


Cl 996 T2002 TS 


Soil 

D, 

o 

D p 

a 


a 

Coarse sand 

690 

1.6 

710 

2 

690 

1.6 

Fine-medium sand 

210 

1.6 

160 

2 

210 

1.6 

Silt 

125 

1.8 

15 

2 

125 

1.8 

Clay 

- 

- 

2 

2 

2 

2 

Salts 

520 

1.5 

- 

2 

520 

1.5 


a Used for the dust flux calculations with the mean mass median diameter 
D p in micrometer. Cl 996 refers to Chatenet et al. [1996], T2002 to Tegen 
et al. [2002], and TS to this study. 
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USGS landuse classes 
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Figure 2. (top) The USGS land use classes (from 1 to 
23) over the domain: (1) urban and built-up land, (2) dry- 
land cropland and pasture, (3) irrigated cropland and pasture, 
(4) mixed dry land/ irrigated cropland and pasture, (5) crop- 
land/grassland mosaic, (6) cropland/woodland mosaic, (7) 
grassland, (8) shrubland, (9) mixed shrubland/grassland, 
(10) savanna, (11) deciduous broadleaf forest, (12) decidu- 
ous needleleaf forest, (13) evergreen broadleaf forest, (14) 
evergreen needleleaf forest, (15) mixed forest, (16) water 
bodies, (17) herbaceous wetland, (18) wooded wetland, (19) 
barren or sparsely vegetated, (20) herbaceous tundra, (21) 
wooded tundra, (22) mixed tundra, and (23) bare ground 
tundra, (bottom) The relative percentage of USGS class 19 
surface (for barren soils). This class is used as a desert mask 
for the dust fluxes calculation when using the USGS data. 


of the five most important land use types in 1° x 1° cells. 
The erodible land use types considered were scrublands (8) 
and barren soil (19). This may appear as a limiting choice, 
but this corresponds to the most important land use classes 
able to produce dust in Northern Africa. This constitutes the 
desert mask as presented in Figure 2. 


2.3. Roughness Length 

[ 19 ] Models use z 0 to calculate surface boundary con- 
ditions for wind by representing the effect on turbulent 
transport of surface elements ranging from vegetation and 
small-scale obstacles to sub-grid orography. Depending on 
the spatial scale, the surface roughness varies widely. In 
meteorological modeling, z 0 is mainly used to estimate the 
friction velocity, a turbulent parameter designed to estimate 
the wind speed near the surface. z 0 in a grid cell will be 
defined such that the modeled wind is realistic compared to 
measurements. In this sense, zo in a model has no “physi- 
cal” and “realistic” meaning: this is a parameter defined to 
integrate correctly the vertical wind. 

[ 20 ] In dust modeling, the roughness length may be used 
to estimate the threshold friction velocity. In meteorological 


and dust models, three roughness lengths at different scales 
may be distinguished as illustrated in Figure 3. 

[ 21 ] 1. At the synoptic scale : This “aerodynamic” rough- 
ness length is used mainly for the atmospheric part of the 
model to calculate the friction velocity to take into account 
the role of subgrid orography, forests, or urbanization when 
slowing the mean wind fields at large scale. Values are often 
look-up tables or just constant over desert surfaces. The val- 
ues generally range from 1 cm to 1 m. This is the case 
of the USGS values presented in Figure 4. This data set is 
not appropriate for the calculation of the threshold friction 
velocity in dust models since it does not represent the rough- 
ness at the scale needed for the emission process over arid 
areas. This data set is used in WRF model (http://wrf- model, 
org/index.php) and is available with the WRF distribution. 

[ 22 ] 2. At the mesoscale : The roughness length accounts 
for small vegetation and obstacles. Values are generally 
lower than the “synoptic” roughness length and range from 
millimeter to centimeter. This is the case of the LISA data 
set: zo is provided after the compilation by Chatenet et al. 
[1996] and Callot et al. [2000]. The values are representative 
of local effects and constitute a more realistic zo data set for 
this study. Another available data set, hereafter called ERS, 
is derived from ERS-1 satellite measurements (with a hori- 
zontal resolution of 0.25° x 0.25°) [ Prigent et al. , 2005]. The 
retrieval approach differs from the in situ based approach 
used for the LISA data set and constitutes an interesting 
alternative for dust emission calculations at a global scale. 

[ 23 ] 3. At the local scale : The smooth roughness length, 
zo, , accounts for mobile erodible sand particles over arid and 
semiarid areas. This length is generally tabulated or calcu- 
lated using the mean mass median diameter of the erodible 
particles, D p , as z 0s = D p / 30, \Marticorena et al, 1997]. 

[ 24 ] In reality, two different types of wind profiles are 
known for fully developed turbulent flow: the wind pro- 
file over aerodynamically rough and over aerodynamically 
smooth surfaces. In the first case, z 0 is a function of the 
roughness elements. At the same time, over arid regions, 
roughness increases apparently due to the suspension of an 
increasing number of particles. At low wind speeds, arid sur- 
faces (also water surfaces) become aerodynamically smooth, 
and z 0 ,' is embedded in the viscous sublayer and scales with 
the kinematic viscosity. 

[ 25 ] In this study, the roughness length is prescribed to 
estimate dust emissions, i.e., we do not take into account vis- 
cous sublayer effects or the apparent increase of roughness 
due to saltating particles. We note that in some advanced 
atmospheric models widely used in the community (e.g., 
NCEP/Eta, NMM5, WRF/ARW, WRF/NMM), viscous sub- 
layers are explicitly taken into account, and the distinction 
between aerodynamically smooth and rough surfaces is thus 
available [Janjic, 1994; Zilitinkevich et al., 2008]. With 
respect to dust emissions, to our knowledge, only the study 
of Nickovic et al. [2001] intends to include viscous sublayer 
effects in a dust emission scheme. 

[ 26 ] Flere, the roughness effect on dust emissions is 
parameterized through a drag partition scheme which 
expresses the efficiency with which drag is partitioned 
between the roughness elements characterized by the aero- 
dynamic roughness length (zo) and the erodible surface 
characterized by a smooth roughness length (z 0v ). The dif- 
ferent values of z 0 (cm) used in this study are presented in 
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Synoptic scale Meso-scale Local scale 

Model Ax=10-1000km Ax=l-10km Ax=l-10()m 

resolution ^ 



Figure 3. The roughness length definitions as a function of its use in local to global meteorological and 
dust transport models. 


Figure 4. The LISA values show an important spatial vari- 
ability and a large range of values (from zq= 1 0 4 to 10 cm). 
The ERS values are spatially variable and close to the LISA 
values. The USGS values are homogeneous over the whole 
domain and are related to the most present land use, barren 
soil. This leads to very high values (1 order of magnitude 
higher than those of LISA and ERS): in this case, a drag par- 
tition scheme is not usable [ MacKinnon et al., 2004], and the 
USGS data set will only be used for the calculation of the 
friction velocity. 


with the constant parameters a„ = 0.0123 and y = 
300 kg nr 2 . The particle density, p p = 2.65 1 0 3 kg m 3 is 
chosen to be representative of quartz grains. 

[ 31 ] The threshold friction velocity, ul, is expressed as 
follows: 


ul(D p ) 


ul s (D p ) x q 

/efK-0, -0s) 


(3) 


with q a soil moisture correction. For the idealized study in 
section 4, this correction is set to q = 1 . 


3. Dust Production Model Design 

[ 27 ] The dust emission fluxes are calculated using the 
parameterization of Marticorena and Bergametti [1995] 
(hereafter denoted MB95) for saltation and the DPM 
proposed by Alfaro and Gomes [2001] (hereafter denoted 
AG01) for sandblasting. Before calculating the fluxes, the 
estimation of the friction velocity is required. All steps 
involved in the calculation are presented in this section. 

3.1. The Friction Velocity 

[ 28 ] The friction velocity, u*, is estimated under neutral 
conditions, as follows: 


Ut 


|U| 


k 

ln(z/zo) 


( 1 ) 


with |(/| the 10 m mean wind speed, k = 0.41, the Kannan 
constant, z the height above ground level where the wind 
speed is estimated by the meteorological model, in this case 
z = 10 m, and zo the roughness length. This roughness length 
is a key point for this study. The same calculation is used 
to estimate dust fluxes in the idealized case and for the 
real case. 

3.2. The Threshold Friction Velocities and the 
Drag Efficiency 

[ 29 ] The saltation flux Fj, is nonzero only if u* > ul(D p ) 
for a given soil particle diameter D p . ul is the threshold 
friction velocity depending on D p , the soil particle diame- 
ter, z 0 , the “aeolian” roughness length, and zo. v , the “smooth” 
roughness length. 

[ 30 ] The “smooth” threshold friction velocity ul s is esti- 
mated following Shao and Lit [2000]: 



40" 

35° 

30" 

25° 

20 ° 

15 ° 

10 ” 



340° 350° 0° 10° 20° 30° 40" 50" 60° 70" 80" 


0.08 

0.07 

0.06 

0.05 

0.04 

0.03 

0.02 

0.01 

0.00 


ERS (OBS-P) 




Figure 4. Comparison of the aerodynamic roughness 
lengths, z 0 (cm), provided with the three surface data sets and 
for the most important relative surface fraction for each cell. 
Horizontal resolution is 1° x 1°. The color scale is the same 
as in Prigent et al. [2005] for LISA and ERS but different 
for USGS, the values being 1 order of magnitude higher. 
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Figure 5. The drag efficiency calculated for a wide range 
of roughness length ad for specific values of the smooth 
roughness length, depending on the mean mass median 
diameter of the particles. 


[ 32 ] f, ff represents the drag efficiency, i.e., a function 
depending on z 0 , the “aeolian” roughness length and z 0s , the 
“smooth” roughness length which is expressed as follows 
[Marticorena et ah, 1997]: 


/eff- 


1-1 


ln(z 0 /z 0s ) 


>(0.35[(0.1/z 0s ) 0 - 8 ]) 


)] 


(4) 


[ 33 ] This function takes into account the non-erodible ele- 
ments dissipating a part of the wind momentum that will not 
be available for saltation. 

[ 34 ] The evolution of this parameter as a function of the 
roughness length values is displayed in Figure 5. The drag 
efficiency is limited to have values between 0 and 1. For 
feg = 1, the threshold friction velocity id is equal to the 
smooth threshold friction velocity ul s : the surface is smooth, 
and there are no non-erodible elements. On the other hand, 
for an increasing roughness length z 0 , non-erodible elements 
are present: / eff tends to zero, increasing the value of id 
and thus limiting the potential erodibility of the considered 
surface. 


3.3. The Wind Speed Distribution 

[ 35 ] Even if the current study uses prespecified wind 
speed over the whole domain in section 4, the goal is to eval- 
uate the mineral dust fluxes under realistic conditions. The 
model grid resolution is 1° x 1°, thus larger than the subgrid 
scale variability of observed winds. 

[ 36 ] If the use of a mean wind speed value may have a 
limited impact on transport, this is not the case for saltation, 
this process being a wind-dependent threshold process. In 
order to taken into account this spatial subgrid scale wind 
variability, a Weibull distribution is implemented, following 
Cakmuretal. [2004] and Pryor etal. [2005]. The probability 
density function is expressed as follows: 


pm) = 

A 



(5) 


where k is a dimensionless shape parameter (in our case 
k = 4, a value commonly used to represent the wind vari- 
ability over large surfaces). A corresponds here to the mean 
modeled wind speed, in meters per second. In our case, this 
mean value corresponds to an hourly value, representative of 
the whole model grid cell. 

3.4. The Horizontal Saltation Fluxes 

[ 37 ] The parameterization of Marticorena and Bergametti 
[1995] estimates saltation fluxes using wind speed values 
and a complete description of the soil characteristics (dry 
size distribution of the loose erodible soil aggregates, soil 
humidity, vegetation cover, presence of non-erodible ele- 
ments such as rocks on the surface). The calculation of the 
vertically integrated saltation flux is based on White’s [1986] 
equation: 


F h (D p ) = K^ui (l - + “*) (6) 

[ 38 ] Following the recommendation of Gomes et ah 
[2003], we use K = 1 in our calculations. The air density 
is also considered as constant (p air = 1.227 kg nr 3 ) for the 
idealized experiments in section 4. 

[ 39 ] Finally, the total saltation flux is obtained by integrat- 
ing F/, over the soil size distribution from D""" = 1 0 6 m to 
£jmax _ 2 1 () 2 m This interval is chosen in order to cover 
the whole range of possible soil sizes. The final saltation flux 
is expressed as follows: 

t- D T x 

F h = / F h (D p )dS K ](D p )dD p (7) 

Jd™ n 

Sxa(D p ) being the relative surface distribution covered 
by particles with a mean mass median diameter D p 
[ Marticorena and Bergametti, 1995]. 

3.5. The Vertical Sandblasting Fluxes 

[ 40 ] In order to estimate size resolved vertical dust fluxes, 
i.e., sandblasting fluxes, the Alfaro and Gomes [2001] DPM 
is used with the numerical optimization described in Menut 
et ah [2005]. The sandblasting flux is computed based on 
the partitioning of the kinetic energy of individual saltat- 
ing aggregates and the cohesion energy of the populations 
of dust particles. This model assumes that dust emitted by 
sandblasting is characterized by three modes whose pro- 
portion depends on the wind friction velocity. From wind 
tunnel measurements performed on two natural soils from 
semi arid regions, Alfaro et ah [1998] consider these three 
modes as independent of the soil types. They described the 
three modes using log-normal distributions with diameters 
rfi = 1.5 x KU 6 m, di = 6.7 x ICE 6 m, and J 3 = 14.2 x 1CT 6 m 
and their associated standard deviation, respectively 
o 1 = 1.7, 02 = 1.6, and 0-3 = 1.5. Based on this model, as 
soil aggregate size or wind speed increases, kinetic energy 
becomes able to release first particles of the coarsest mode 
that are associated with the lowest cohesion energy, then par- 
ticles from the intermediate population, and finally the finest 
particles. It also implies that for a specific wind speed and 
soil size distribution, the dust flux may be zero even if the 
saltation process occurs. In order to apportion the available 
kinetic energy between the three modes, a constant cohesion 
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Figure 6. Main African dust sources, (top) The aerosol 
optical depth at 555 nm measured by MISR and averaged 
from February 2000 to December 2010, with a 0.5° x 0.5° 
resolution, (bottom) The Goddard Chemistry Aerosol 
Radiation and Transport (GOCART) preferential sources, 
after Ginoux et al. [2001]. 


energy e, is associated to each mode values. The numeri- 
cal values of e, were determined by adjusting the predicted 
aerosols size distribution to those measured in wind tun- 
nel under different wind conditions, using an iterative least 
square routine. In this study, the values recommended by 
Alfaro and Gomes [2001] are used: e\ = 0.376, e 2 = 0.366, 
and e 3 = 0.346 kg m 2 s 2 . 

[ 41 ] The kinetic energy is expressed as a function of the 
soil particle diameter after Alfaro et al. [1997] and Shao and 
Lu [2000]: 


100tt 


Pp~ 


-Dliu.f 


( 8 ) 


[ 42 ] It is compared to the cohesion energy of the three 
aerosol modes in order to compute the proportion pt(D p ) of 
these three modes in the total dust size distribution. By inte- 
grating equation (8) over the three aerosols modes, the total 
sandblasting flux may be written as follows: 


Fv,m,i(Dp ) 


J D™' n 


~Pp\ 


, Ri(D„,;)4. 


d F h (D p f 


(9) 


where i= I ,/Y c iass is the number of intervals discretizing the 
soil size distribution in the range [/J"’ in : /.T m,x ] and d m i the 
mean mass diameter [Menut et al., 2005]. 


4. Idealized Experiments 
4.1. Experimental Setup 

[ 43 ] The sensitivity of the modeled mineral dust fluxes 
to different surface and soil descriptions is estimated 
over the simulation domain covering Northern Africa and 
Middle East. In this section, we consider spatially constant 
wind speeds for the analysis. The horizontal resolution was 
chosen to be the same for all calculations, 1° x 1°. For 


each of these configurations, the Alfaro and Gomes [2001] 
DPM is used, after its numerical optimization described in 
[Menut et al. Menut et al. [2005]]. We define three basic 
experiments: 

[44] 1 . NOfg : without the drag efficiency calculation and 
thus with / ef f forced to 1 . 

[45] 2. Feff. with the drag efficiency as described in 
equation (4). 

[ 46 ] 3. GOpref with a preferential source function. This 
approach was proposed partly due to the historical lack of 
roughness length [ Ginoux et al., 2001] and has been widely 
used to specify the most erodible areas in models [Cakmur 
et al., 2006; Perez et al., 2011]. In Figure 6 (bottom), the 
values represent a multiplicative factor, a, with a horizon- 
tal resolution of 1° x 1°, ranging from 0 to 1 and built to 
enhance the representation of the most erodible surfaces. In 
this case, the vertical dust fluxes are calculated using/ eff = 1 
and F i cioprci et, x F v . 

[47] In addition to each of the previous experiments, we 
define three different uses of the roughness length and the 
soil types, as described in Figure 7. 

[ 48 ] 1. LISA: represents the flux with the use of the LISA 
data set for the soil texture and the roughness length in the 
drag partition scheme. 

[49] 2. ERS-LISA: represents the flux using the LISA-soil 
texture and the ERS satellite zo in the drag partition scheme. 

[50] 3 . ERS-SGF : represents the flux using the STATSGO- 
FAO for the soil texture and the ERS satellite z 0 in the drag 
partition scheme. 

[51] Note that the soil texture affects the threshold fric- 
tion velocity through z 0 , which depends on the coarser soil 
particle population. It also affects the saltation flux through 
the relative surface area of the soil particles and the kinetic 
energy in the vertical flux though the soil particle diameter. 

[ 52 ] The previous combinations lead to the calculation 
of nine different mineral dust fluxes. These nine configura- 
tions make possible to evaluate the relative sensitivities of 
the fluxes to the soil data sets and to the way of prescibing 
roughness over erodible surfaces. Thecomparisons between 
LISA and ERS-LISA are used to estimate the impact of the 
roughness length on the flux calculations. The comparisons 
between ERS-LISA and ERS-SGF are used to quantify the 
impact of the soil texture database on the flux calculations. 



Figure 7. General flowchart of the model configurations 
used in this study. See details in section 4. The “Feff’ and 
“GOpref’ are optional and may be activated or not during 
the calculation. 
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Table 3. Names and Coordinates of the Main Dust Sources Areas 3 


Area Name 

Longitude 

Latitude 

Whole domain 

-16°W -+ +77°E 

+13°N -> +36°N 

Bodele 

+I1°W — + +19°W 

+15°N— * +19°N 

Mali 

-10°W -+ -2°W; 

+17°N — + +21°N 

Chotts 

+5°W +10°W 

+29°N ->• +35°N 

Lybian Desert 

+17°W —*■ +25°W 

+23 °N -+ +30°N 

Arabian Desert 

+45°W -+ +54°W 

+18°N — + +24°N 


“Defined for the results presentation and displayed in Figure 6. 


[53] The main areas of analysis are displayed in Figure 6 
(top) where the map represents the aerosol optical depth 
(AOD) at 555 nm measured by Multiangle Imaging Spec- 
troradiometer (MISR) averaged from February 2000 to 
December 2010 at 0.5° x 0.5° resolution. The MISR AOD 
represents both sources and dust load in the atmosphere after 
transport. However, it can be considered that the strongest 
regional AOD roughly coincides with the main dust sources. 
Among the sources, we highlight the well-known hot spots 
in the Bodele and eastern Niger, the Mali/Mauritania bor- 
der, the Chotts, the Lybian desert, and the Arabian desert. In 
Figure 6 (top), the rectangles mark the areas for which the 
results are presented. Results are also presented for a larger 
domain encompassing the whole Saharan-Sahel region. This 
domain corresponds to the bounds within which the LISA 
data are available. Names and coordinates of the areas are 
given in Table 3. 

4.2. Results and Discussion 
4.2.1. Horizontal Maps 

[54] Maps of surface fluxes for the nine configurations are 
displayed in Figure 8 for a spatially constant wind speed of 
12 m/s. Note that all results are presented as a function of 
wind speed in contrast to other studies using friction velocity 
[e.g., Marticorena and Bergametti, 1995; Darmenova et al., 
2009]. This choice is required for the present study since 
we analyze the impact of the roughness length which is also 
used in the calculation of the friction velocity. 

[55] The NOfg experiment (Figure 8, left column) cal- 
culates emission fluxes which are considerably higher than 
those derived from the Feff and GOpref experiments. Also, 
the spatial variability is rather low over the whole domain 
displaying a poor correspondence with the dust hot spots 
detected by MISR. The worst configuration for this exper- 
iment is certainly LISA, which presents unrealistic fluxes 
over the Sahel and the northern part of the domain. When 
using the ERS roughness length, the results are more realistic 
with both LISA and STATSGO-FAO soils. 

[56] For the Feff experiment (Figure 8, middle column), 
the LISA configuration represents fairly well the spatial vari- 
ability of dust emission. ERS-LISA retains the major sources 
but appears to be more restrictive than LISA. ERS-SGF 
shows a close spatial variability to ERS-LISA due to the use 
of the same roughness length. However, the flux magnitudes 
are higher due to the different soil texture, which is smoother 
for the STATSGO-FAO soil than for the LISA soil. 

[57] For the GOpref experiment (Figure 8, right col- 
umn), the three configurations exhibit a rather close spatial 
variability showing that the fluxes in this case are mostly 
constrained by the a factor. In contrast to the Feff exper- 
iments, the preferential source enhances the fluxes in the 


border of Algeria with Mauritania and Mali and omits to a 
large extent the Mali/Mauritania border source. 

[58] The results reflect that the use of a drag partition 
scheme strongly constrains the spatial distribution of fluxes 
and better matches the satellite dust distribution. The fluxes 
are also lower by a factor of 2 to 3 when using the drag 
partition scheme. Also, for a given soil data set, the fluxes 
estimated using the two roughness length data sets show sig- 
nificant differences in the spatial distribution, while for a 
given roughness length data set, the fluxes estimated using 
the two soil texture data sets show significant differences in 
the intensity of the fluxes. These points will be discussed in 
more detail in the next section. 

4.2.2. Distributions of Fluxes 

[59] Figure 9 shows a probability distribution of dust 
fluxes for the NOfg, Feff, and GOpref experiments. The flux 
values are represented with a step of 0.5 1 0 8 g cnr 2 s 1 . The 
occurrence refers to the fraction of grid cells where each dust 
flux value is estimated. 

[60] All configurations show a strong positive skew 
(towards lower flux values). However, for all experiments, 
we notice that LISA displays significantly different distri- 
butions than ERS-LISA and ERS-SGF. With the LISA zq, 
fluxes are less positively skewed and exhibit pronounced 
peaks in contrast to the smoother distributions observed in 
ERS-LISA and ERS-SGF. 

[61] The NOfg experiment is less positively skewed than 
Feff and GOpref experiments. It clearly appears that Feff and 
GOpref act on the low flux values by decreasing the mode of 
the distributions from 10-15 x 10 8 to 1-3 x 10 s g cm 2 s ! . 

[62] The distributions for the Feff experiment reveal the 
impact of the soil and the roughness data sets on the intensity 
of the fluxes, with the roughness length data set having a 
stronger effect than the soil data set. 

[63] The more comparable distributions are obtained with 
the GOpref experiment. The weight of the preferential 
source constrains the distributions to similar values between 
0 and 5 x 10~ 8 g cnr 2 s 1 with the ERS zo and between 0 and 
15x1 0 8 g cnr 2 s 1 with the LISA z 0 . 

4.2.3. Relative Impact of z 0 Distributions 

[64] Figure 10 presents distributions of aeolian rough- 
ness length, zo, over the whole African domain for the two 
data sets, LISA and ERS. Since the LISA data set is based 
on geomorphological considerations and measurement data, 
z 0 values are prescribed for some specific land types, and 
thus, only “discrete” values are available. On the other hand, 
the ERS zo values are estimated from satellite measure- 
ments, exhibiting a continuous spectrum of values and a 
maximum zo of 0.08 cm. It is obvious that, to describe phys- 
ical characteristics, a continuous spectrum is more realistic 
than some few discrete values, and in the case of min- 
eral dust flux calculations, this leads to large differences 
since the parameterizations used are based on a threshold 
approach. 

[65] To better quantify this “discrete values” effect, 
Figure 11 presents a probability distribution of the drag 
efficiency, / eff , for LISA, ERS-LISA, and ERS-SGF at the 
prespecified wind speed of 12 m/s. Occurrence is calculated 
with a step of 0.0 1 . This representation explains the main dif- 
ferences diagnosed in the flux calculations. Since the LISA 
z 0 data set used only a few different values of z 0 , the corre- 
sponding / eff also exhibit a few discrete values. On the other 
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Figure 8. Mineral dust vertical fluxes (10 8 g cur 2 s ') calculated with a prespecified wind speed of 
12 m/s with the model configurations LISA, ERS-LISA, and ERS-SGF, (left) with no drag efficiency or 
preferential source use, called NOfg, (middle) using the drag efficiency f e s scheme, and (right) using the 
GOCART preferential sources, noted GOpref. Note that the color scale is exponential. 


hand, when f eS is estimated with the ERS z 0 , the / e ff is more 
continuous and more physically realistic. Also, there is a 
minimum drag efficiency with the ERS zq of about 0.3 which 
explains the more restrictive regions of emission displayed 
in Figure 8. 

4.2.4. Fluxes Over Specific Areas 

[66] For all configurations, the dust production model was 
run with spatially constant wind speed values, from 0 to 
24 m/s with a step of 1 m/s. Results were spatially integrated 
over specific source regions, and the results are displayed in 
Figure 12. An integration of fluxes over the whole domain 
simulation is referred to as “ALL.” 

[67] The fluxes with the configuration LISA tend to be 
higher than those obtained with the ERS z (l (with both ERS- 
LISA and ERS-SGF). This is observed for the whole domain 
(“ALL”) integration and for some specific regions, even if 
this tendency is lower over Bodele and the Lybian and Ara- 
bian deserts. Another general tendency is that the fluxes 
without drag efficiency (NOfg, GOpref), are significant at 
very low wind speeds (2-6 m/s). In arid regions, modeled 
emissions fluxes are often close to zero below 10 m wind 
speeds of 7-8 m/s [Menut et al, 2005]. This effect is clearly 
emphasized with the LISA configuration and is certainly 
unrealistic. The same applies for ERS-LISA and ERS-SGF 
but to a lesser extent. Globally, the omission of a drag 
partition scheme gives unrealistic emission fluxes for wind 
speeds between 4 and 8 m/s. For wind speed values up to 
8 m/s, the configurations NOfg and GOpref produce higher 
fluxes than Feff for all regions. 

5. Real Case Study 

[68] In order to evaluate model configurations with real- 
world data, we simulated the period 1 March to 31 July 
2011 using the WRF and CHIMERE models. The modeled 
dust concentration fields were used to estimate the Aerosol 


Optical Depth (AOD) which was compared to Aerosol 
Robotic Network (AERONET) Sun photometer measure- 
ments. In this section, we first describe the dust transport 
model and the observations used. Results are presented 
under the form of time series and AOD maps and are 
discussed using skill scores. 

5.1. Model Setup and Observations 

5.1.1. The WRF and CHIMERE Models 

[69] Two models are used in this study: the meteorolog- 
ical fields are calculated using the Weather Research and 
Forcasting (WRF) model [Skamarock et al., 2007], and the 
mineral dust transport is calculated with the CHIMERE 
chemistry-transport model in its mineral dust configuration 
[Menut et al, 2009]. 

[ 70 ] The meteorological analysis data of NCEP [Kalnay 
et al., 1996] are provided on a regular 1.125° x 1.125° grid. 
These fields are used as boundary conditions and nudging 
for WRF in its 3.2.1 version. The model is used in its non- 
hydrostatic configuration. The horizontal grid of 180 x 102 
grid points in longitude and latitude is used with a Lambert 
projection to ensure a regular grid spacing of Ax = A y = 
100 km. The vertical grid covers 32 levels from the surface 
to 50 hPa, and the integration time step is 4 min. For the 
microphysics, the WRF Single Moment-5 class scheme is 
used allowing for mixed phase processes and supercooled 
water [Hong et al, 2004]. The radiation scheme is RRTMG 
scheme | M lower et al, 1997]. The surface physics is cal- 
culated using the Noah Land Surface Model scheme with 
four soil temperature and moisture layers [Chen and Dudhia, 
2001]. The planetary boundary layer physics is processed 
using the Yonsei University scheme [Hong et al, 2006], and 
the cumulus parameterization uses the ensemble scheme of 
Grell and Devenyi [2002]. 

[ 71 ] The CHIMERE chemistry-transport model may 
be used for urban to regional air pollution (including 
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Figure 10. Distribution of eolian roughness lengths, z 0 (in 
cm), over the whole African domain for the LISA and ERS 
data sets. 

is the one of Troen and Mahrt [1986]. Dry and wet 
deposition are treated as described in Wesely [1989] and 
Loosmore and Cederwall [2004]. The horizontal resolution 
is 1° x 1° with 20 vertical levels (hourly dependent on 
the pressure) from the surface to 200 hPa. The calcula- 
tion time step is 10 min, and the results are stored every 
hour for analysis. Model characteristics are explained in 
Menut et al. [2007]. In order to be realistic for this simu- 
lation, the configuration Feff is used for the three soil and 
surface data sets. 

5.1.2. AERONET Optical Thickness Measurements 

[ 72 ] The AERONET stations [Holben et al, 2001] used 
are listed in Table 4. Figure 13 shows the stations’ loca- 
tions. Data are extracted from the NASA website in their 
1.5 version that includes cloud screening. The retained sta- 
tions were chosen to represent both sites close to important 
sources (Banizoumbou, Cinzana, Ilorin, Dakar) and sites 
more representative of midrange dust transport (Capo Verde, 
Izana, Blida, and Lampedusa). 

5.2. Results 

[ 73 ] Figure 13 presents horizontal maps of AOD. Hourly 
values are averaged over 5 days and only for hours 9:00, 


Figure 9. Distribution of vertical dust fluxes calculated 
with several model configurations and for a prespecified and 
constant wind speed of 12 m/s. The calculations are done 
with a step of 0.1 x 1 0 x g cnr 2 s 1 for the dust flux. The 
occurrence represents the number of grid cells where each 
dust flux value was estimated. 


anthropogenic and biogenic emissions [Bessagnet et al, 
2010]), for analysis and forecast [Menut and Bessagnet, 
2010], and for long-range transport [Menut et al., 2009]. 
In this study, the model is used with only the mineral dust 
emissions and transport. The domain covers the area (-56° 
< longitude < +38°) and (+2° < latitude < +54°). It takes 
into account the processes that drive the main phases of the 
dust cycle: saltation and sandblasting for emissions calcula- 
tions, horizontal advection, and vertical transport, including 
advection and mixing. The boundary layer parameterization 



Figure 11. Distribution of the drag efficiency /di calculated 
for the three model configurations and a prespecified wind 
speed of 12 m/s. 
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Region: ALL + NOfg Region: ALL + Feff Region: ALL + GOpref 








Figure 12. Mineral dust fluxes (1(T 8 g cur 2 s 1 ), accumulated over specific source regions and for the 
three model configurations LISA, ERS-LISA, and ERS-SGF. (left) With no drag efficiency or preferential 
source used, (middle) using the drag partition scheme, and (right) with the GOCART preferential source. 


10:00, and 11:00 to produce one single map comparable to 
MODIS Terra satellite data. 

[ 74 ] Two different periods are presented: 30 April to 4 
May 2011 and 6-11 June 2011. The top panel corresponds 
to the MODIS satellite aerosol optical at 550 nm. The 
modeled AOD was also calculated at 550 nm for this com- 
parison to MODIS. The other panels correspond (from top 
to bottom) to the LISA, ERS-LISA, and ERS-SGF model 
configurations. For the two periods, and even if no satellite 
data are available over the middle of Africa, the observed 
AOD highest values are well localized with the model. For 
the first period (30 April to 4 May), the main active sources 
are the Lybian desert, Bodele, and the Chotts region. The 
comparison between LISA and ERS-LISA shows similar 
spatial patterns with the dust content being higher towards 
the Mediterranean. As shown in the previous section, the 
ERS zo provides more limited emission areas for a given 
wind speed, and the maxima are lower. The relative patterns 


between north and south significantly differ with ERS-SGF 
compared to the LISA configurations. In this case, the 
Bodele region shows higher AOD than the Mediterranean 
region. 

[ 75 ] For the stations listed in Table 4, AOD time series 
at 440 pm are compared with AERONET data, Figure 14. 
At Banizoumbou (Nigeria), Zinder (Nigeria), and Cinzana 


Table 4. Names and Coordinates of the AERONET Surface 
Stations Displayed in Figure 13 


Station 

Longitude 

Latitude 

Banizoumbou 

+2.66°E 

+13.54°N 

Capo Verde 

-22.93°W 

+16.73°N 

Cinzana 

-5.93°W 

+13.28°N 

Dakar 

— 16.95°W 

+14.39°N 

Lampedusa 

+12.63°E 

+35.52°N 

Zinder 

+8.98°E 

+13.75°N 


6515 


MENUT ET AL.: IMPACT OF ROUGHNESS AND SOIL TEXTURE ON MINERAL DUST 



MODIS Terra April, 30 to May, 4 2011 
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LISA April, 30 to May, 4 2011 
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Figure 13. Mean averaged values of AOD during periods of 5 days, starting with 30 April 2011 and 6 
June 2011. The MODIS Terra optical depth at A = 550 nm are extracted from the NASA Giovanni tool. 
The LISA, ERS-LISA, and ERS-SGF are modeled aerosol optical depth (A = 550 nm) with WRF and 
CHIMERE. 


(Mali), which are close to mineral dust sources in the 
south, the AOD values are higher than over the other sites, 
with measured peaks reaching values up to 2. While the 
LISA and ERS-SGF configurations reproduce the observed 
variability between 0 and 2, the ERS-LISA configura- 
tion AOD remains below 0.5. This remark is valid for 
all sites. 

[76] Dakar is located at the West African coast and 
Sal/Cape Verde about 600 km off the West African coast. 


For these two locations, the LISA configuration estimates 
higher concentrations than the two other configurations and 
the observed AOD. The ERS-SGF configuration is the one 
that better captures the AOD at these two stations. At the 
Lampedusa station (Italy, 150 km off the Tunesian coast), the 
measured AOD values are rather low (0.5 in average) and 
representative of long-range transport of dust. Here, the three 
configurations show a strong variability, and it is difficult to 
qualitatively conclude on the most suitable configuration. 
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Figure 14. Time series of aerosol optical depth (A = 440 nm), measured with the sumphotometers of 
the AERONET network and modeled with CHIMERE and the several dust emissions data sets. 


[ 77 ] Statistical scores are calculated over the period of 1 
March to 1 July 2011. Table 5 presents these results with the 
following statistical indicators for each station: 

[ 78 ] 1 . Correlation: 


Cor _ E^W-M)(o,-o) 

[ 79 ] 2. Mean normalized bias error (MNBE): 


( 10 ) 


MNBE = - V ( M ‘ x 100 (11) 

JVtfV O t J 


[so] 3. Root mean square error (RMSE): 


RMSE 


i V(AC-0,) 2 


\| 


( 12 ) 


where M and O represent the modeled and mesured AOD, 
respectively, with an hourly frequency (corresponding to 


the model output). M and O represent the time-averaged 
AOD for the model and the measurements, respectively. 
The N values correspond to the number of values for which 
AERONET data are available, which differs from one sta- 
tion to another. The best values of correlation and MNBE 
are in bold for each station to show the most accurate model 
configuration. 

[si] The mean bias (MNBE) is expressed in percent and 
ranges from -73% to +175%. This indicator has to be care- 
fully analyzed because the bias is directly dependent on the 
model dust size-distribution discretization (the number of 
bins and their values) and on the distribution of emissions in 
the model modes. It means that a correlation is more model 
dependent than a real indicator of a good (or not) modeling. 
But this may globally express the range between a model and 
observations. For the three model configurations, the biases 
are quite comparative and mainly show that all model con- 
figurations are able to give the correct order of magnitude 
for the AOD. 


Table 5. Statistical Scores for the AOD Modeling 3 


Site 

^val 

O 


M 



RMSE 



Cor 



MNBE 


LI 

EL 

ES 

LI 

EL 

ES 

LI 

EL 

ES 

LI 

EL 

ES 

Banizoumbou 

1068 

0.66 

0.47 

0.17 

0.74 

0.41 

0.57 

0.56 

0.18 

0.49 

0.50 

-15.73 

-70.99 

14.32 

Zinder 

826 

0.68 

0.46 

0.20 

0.84 

0.43 

0.60 

0.65 

0.36 

0.44 

0.37 

17.71 

-66.11 

39.64 

Cinzana 

583 

0.76 

0.50 

0.18 

0.65 

0.58 

0.74 

0.64 

0.21 

0.32 

0.27 

-19.07 

-73.91 

-6.30 

Capo Verde 

450 

0.26 

0.47 

0.08 

0.19 

0.38 

0.20 

0.14 

0.71 

0.74 

0.68 

67.81 

-70.59 

-27.05 

Dakar 

1132 

0.54 

1.35 

0.18 

0.38 

1.11 

0.44 

0.34 

0.37 

0.45 

0.47 

175.49 

-65.47 

-28.25 

Lampedusa 

671 

0.24 

0.27 

0.05 

0.08 

0.29 

0.23 

0.20 

0.74 

0.68 

0.74 

-3.94 

-81.45 

-72.46 


“Compared to the AERONET measurements (for X = 440 nm). The letters correspond to “LI” for LISA, “EL” for ERS-LISA, and “ES” for ERS-SGF. 


6517 


MENUT ET AL.: IMPACT OF ROUGHNESS AND SOIL TEXTURE ON MINERAL DUST 


[ 82 ] The correlations, “Cor” in Table 5, range from 0.18 
to 0.74, depending on the site and the model configuration. 
This indicator is much appropriate to the discussion of this 
study, because this is the direct expression of the ability of 
different modeled sources to be well represented with a real- 
istic meteorology and after transport in the atmosphere. For 
the LISA configuration, the correlation ranges from 0. 1 8 to 
0.74, for ERS-LISA from 0.32 to 0.74, and for ERS-SGF 
from 0.27 to 0.74. The high values are not estimated for 
the same sites, but it is noticeable that the best correlations 
are calculated for Lampedusa and Capo Verde sites, after 
transport. The correlations shows that the three model con- 
figurations are able to produce realistic transport and AOD. 
This is already known for the LISA data set being used in 
many dust models. With these scores, we show that this is 
also possible to have good model results with a new data set 
such as ERS-SGF. 

[83] The RMSE expresses the model variability compared 
to observations. It varies from 0.29 to 1.11 with LISA, 
from 0.2 to 0.74 for ERS-LISA, and from 0.2 to 0.65 for 
ERS-SGF, showing that the model variabilities are compa- 
rable between the data sets, the ERS-SGF having the lowest 
values. 

6. Summary and Conclusions 

[84] In this paper, mineral dust emission fluxes modeled 
over Northern Africa were compared when using different 
roughness and soil texture data sets. For all comparisons, we 
used the dust production model of Alfaro and Gomes [2001]. 

[85] Numerous dust modeling studies exist for Northern 
Africa and the Sahel regions. One of the most known and 
widely used soil (texture) and surface (z 0 , the roughness 
length) data sets is the one proposed by LISA and described 
in Marticorena and Bergametti [1995], with 1° x 1° hori- 
zontal resolution. Data were estimated from field campaigns 
and geomorphological analyses. 

[86] We analyzed the sensitivity of an advanced emis- 
sion scheme in Northern Africa by comparing the LISA data 
set to global estimates of soil texture (STATSGO-FAO) and 
roughness length (USGS and ERS). USGS values for rough- 
ness length depend on land use types, and although they 
have been widely used in atmospheric modeling, their appli- 
cation is limited for dust modeling purposes. ERS are more 
recent estimates of roughness length from satellite data and 
are more appropriate for dust modeling. 

[87] Three surface and soil data sets are thus defined, fol- 
lowing Figure 7: the original LISA data set, the LISA data 
set but with the ERS roughness length (ERS-LISA), and the 
USGS and STATSGO surface and soil data sets with the ERS 
roughness length (ERS-SGS). I addition, we also analyzed 
the impact of the drag partition efficiency process in the 
dust emission scheme and compared its behavior within the 
model to a widely used preferential source approach [Ginoux 
et al., 2004]. 

[88] First, and in order to quantify the differences between 
modeled emissions fluxes that are only due to the surface 
and soil characteristics, we defined idealized experiments 
by using a prespecified constant wind speed over the whole 
domain. Having no direct emission flux measurements, com- 
parisons were done among the configurations. The use of 
the drag partition scheme is able to fairly well constrain the 


spatial distribution of the sources with both LISA and ERS 
roughness lengths. The LISA dust fluxes are, on average, 
higher than the other model configurations. This is mainly 
due to the distribution of z 0 values: being tabulated, val- 
ues are more discrete and thus less variable and realistic 
than those provided by ERS from satellite data. Concerning 
the soil texture, the STATSGO-FAO soil data set is real- 
istic enough to make possible fluxes calculations over the 
whole domain. It was also shown that without a drag parti- 
tion scheme, fluxes are unrealistic. The use of a preferential 
sources approach replacing the drag partition scheme is able 
to match the spatial distribution of important sources areas. 
However, this approach does not allow constraining real- 
istically the threshold friction velocity with the emission 
scheme used, since fluxes are too high for low wind speeds. 

[89] Second, the three soil and surface data set configura- 
tions were used to model the period of March to July 2011 
over a domain encompassing Northern Africa and south- 
ern Europe. The models used were WRF for meteorology 
and CHIMERE for the dust transport. The simulations were 
compared to AOD measured at several AERONET stations 
in the region and to MODIS satellite data. The evalua- 
tion showed that the new configuration ERS-SGF is able to 
accurately reproduce the main dust sources and the AOD 
variability in the region. 

[ 90 ] The next steps for these studies will be to compare 
modeled dust fluxes calculated with the ERS data used in 
this study (0.25° x 0.25° [Prigent et al., 2005]) to the recent 
data set presented in Prigent et al. [2012] (6 km x 6 km) and 
thus to explore the impact of the model horizontal resolu- 
tion on emission dust fluxes magnitude and size distribution. 
Another perspective is how to use this global data set to 
extend the dust emissions domain to estimate the impact 
of local erosion on the budget of surface particulate matter 
in Europe, following the study of Bessagnet et al. [2008] 
about the dust source in Ukraine. Finally, the ability to 
model dust emissions fluxes at the same time in Africa and 
Europe would improve the pre visibility of particles and thus 
improve the short-term forecast of air quality, as discussed 
in Mennt and Bessagnet [2010]. 
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